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A rapidly developing field of spintronics is 
based on the premise that substituting charge 
with spin as a carrier of information can lead 
to new devices with lower power consump- 
tion, non- volatility and high operational speediiS. 
Despite efficient magnetization detection^i^i^, 
magnetization manipulation is primarily per- 
formed by current-generated local magnetic fields 
and is very inefficient. Here we report a 
novel non-volatile hybrid multiferroic memory 
cell with electrostatic control of magnetization 
based on strain-coupled GaMnAs ferromagnetic 
semiconductor^ and a piezoelectric material. We 
use the crystalline anisotropy of GaMnAs to 
store information in the orientation of the mag- 
netization along one of the two easy axes, 
which is monitored via transverse anisotropic 
magnetoresistance^. The magnetization orienta- 
tion is switched by applying voltage to the piezo- 
electric material and tuning magnetic anisotropy 
of GaMnAs via the resulting stress field. 

In magnetic memories the information is stored in the 
orientation of magnetization. The major weakness of 
current MRAM implementations hes in the inherently 
non-local character of magnetic fields used to flip ferro- 
magnetic domains during the write operation. Neither 
thermal switching^ nor current induced switching^ can 
solve the problem completely. 

An attractive alternative to conventional ferromagnets 
are multiferroic materialsiS, where ferromagnetic and fer- 
roelectric properties coexist and magnetization can be 
controlled electrostatically via magnetoelectric coupling. 
In single-phase and mixed-phase multiferroics ferromag- 
netic material must be insulating in order to avoid short- 
circuits. Alternatively, magnetoelectric coupling can be 
introduced between ferromagnetic and ferroelectric ma- 
terials indirectly via strain, and in this case the ferromag- 
netic material can be conducting. A conceptual memory 
device made of ferromagnetic and piezoelectric materials 
has been proposed^, and permeability changes in mag- 
netostrictive filmai^, changes in coercive field^'^ and re- 
orientation of the easy axis from in-plane to out-of-plane 
in Pd/CoPd/Pd trilayers^'' have been recently demon- 
strated. 

Some ferromagnetic materials have a complex 
anisotropic magnetocrystalline energy surface and sev- 
eral easy axes of magnetization. If the energy barrier 
between adjacent easy-axis states can be controlled by 
strain, then a multiferroic non-volatile multi-state mem- 




FIG. 1: Top: sketch of the Hall bar with relative orientation 
of electrical current /, magnetic field H and magnetization M. 
Bottom: an AFM image of a 2/jm-wide Hall bar is schemati- 
cally shown attached to the piezoelectric. Strain is applied to 
the sample along [100] and [010] directions, red dashed lines 
on the sketch. 



ory device can be realized. Magnetic anisotropy of dilute 
magnetic semiconductor GaMnAs films is largely con- 
trolled by epitaxial strai n^^'^^ , with compressive (ten- 
sile) strain inducing in-plane (out-of-plane) orientation 
of magnetization. In GaMnAs compressively strained 
to (001) GaAs there are two equivalent easy axes of in- 
plane magnetization: along [100] and along [010] crystal- 
lographic directions. Lithographically-induced unidirec- 
tional lateral relaxation can be used to select the easy 
In addition to the large in-plane crystalline 
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anisotropy there is an uniaxial anisotropy between [110] 
and [110] directions which is probably due to the un- 
derlying anisotropy of the reconstructed GaAs surface^?.. 
We use this magnetocrystalline anisotropy in combina- 
tion with the magnetostrictive effect to demonstrate a 
bi-stable memory device. 
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FIG. 2: Transverse AMR as a function of magnetic field angle 
ifH for unstrained (a) and strained (b,c) samples. Arrows 
indicate the field sweep directions. Magnetic field is _ff = 50 
mT and T = 35 K. Voltage applied to the piezoelectric Vpzt ~ 
for the curves (b) and Vp^t = ±200 V and for (c). 



Molecular beam epitaxy at 265° C was employed to 
grow 15-nm thick epilayers of Gao.92Mno.08As on semi- 
insulating (001) GaAs substrates. The wafers were sub- 
sequently annealed for 1 hour at 280° C in a nitrogen 
atmosphere. Annealing increases the Curie temperature 
of the ferromagnetic film to ~ 80 K and reduces 
the cubic anisotropy. The GaMnAs layer was patterned 
into 2 /im-wide Hall bars oriented along the [110] axis 
by combination of e-beam lithography and wet etching, 
see Fig. [T] After lithography 3 mm x 3 mm samples 
were mechanically thinned to ~ 100 /im and attached to 
a multilayer PZT (piezoelectric Icad-zirconium-titanate 
ceramic)^ with epoxy, aligning the [010] axis with the 
axis of polarization of the PZT. Application of positive 
(negative) voltage VpzT across the piezoelectric intro- 
duces tensile (compressive) strain in the sample along the 
[010] direction, and strain with the opposite sign along 
the [100] direction proportional to the piezoelectric strain 
coefficients d^s ~ —2^31. Both coefficients decrease by a 
factor of 15 between room temperature and 4.2 K. The in- 
duced strain e ~ AL/L for both [010] and [100] directions 
was monitored with a biaxial strain gauge glued to the 
bottom of the piezoelectric. Strain is proportional to the 
change of the gauge resistance, and was measured in the 
Wheatstone bridge configuration: Ae — £[010] ~ £[100] = 
(AL/L)[oio] - (AL/L)[ioo] = "(^[oio] -R[wo])/R, where 
a is the gauge sensitivity coefficient and R is the resis- 
tance of the unstrained gauge. It has been shown before^^ 
that the strain gradient across the piezoelectric and the 
sample is negligible: i.e., gauges glued on top of the sam- 
ple and on the opposite side of the piezoelectric measure 



FIG. 3: Transverse AMR (TrAMR) for the strained Hall bar 
as a function of uniaxial strain measured by a strain gauge. 
Static magnetic field H — 50 mT is applied at ipn = 62°, 
the data is taken at T = 35 K. On the top axis approxi- 
mate voltage on the piezoelectric is indicated. Orientation of 
magnetization relative to the current flow is shown schemat- 
ically for the TrAMR> and TrAMR< states. In the inset 
TrAMR is plotted for J? = 50 mT (red), 70 mT (green) and 
100 mT (blue) measured at T = 25 K. 

similar strain. 

Direction of the in-plane magnetization M is measured 
via transverse anisotropic magnetoresistance (TrAMR), 
also known as giant planar Hall effect^- 

TrAMR^ {pi^ - pj_)sinipm COS ipm, (1) 

where pn and p± are the resistivities for magnetization 
oriented parallel and perpendicular to the current. The 
sign and magnitude of TrAMR depend on the angle 
tpm between magnetization M and current /||[110], see 
schematic in Fig. [TJ TrAMR reaches minimum (max- 
imum) when M||[010] (M||[100]). Longitudinal AMR 
is oc cos^ (fim and is not sensitive to the magnetization 
switching between 45° and 135°. 

In Fig. [2] TrAMR is plotted for the strained and un- 
strained Hall bars as magnetic field of constant magni- 
tude H ^ 50 mT is rotated in the plane of the sample. 
For the unstrained sample (not attached to the piezo- 
electric) there are four extrema for the field angle ipn 
around 45°, 135°, 225° and 315° with four switchings 
of magnetization by 90° near [110] and [110] directions, 
which rcfiects the underlying crystalline anisotropy of 
GaMnAs. For the sample attached to the piezoelectric 
there is a gradual change in the angle of magnetization 
with only two abrupt 90° switchings of the direction of 
magnetization per field rotation, which indicates strong 
uniaxial anisotropy due to a highly anisotropic thermal 
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FIG. 4: Experimentally measured (left panel) and modelled 
(right panel) transverse AMR for the strained sample at dif- 
ferent voltages on the piezoelectric (for comparison data for 
an unstrained sample are shown in the inset). All data were 
taken at T = 35 K and H = 50 mT. In the model we used 
experimentally measured strain Ae, other model parameters 
are discussed in the text. 



FIG. 5: a) Polar plot of magnetostatic energy E/M as a 
function of magnetization angle (prn for H = (black) , 50 mT 
(blue) and 100 mT (red) for ipn = 62° and Vp^t = (Aeif^ = 
16 mT). b-d) Angular dependence of E/M for AeK^ = 13-19 
mT (black to magenta) ior H = (e), H ^ 50 mT (b,d) and 
H — 100 mT (c) . Blue line and arrow mark ipH = 62° , dashed 
red lines indicate two stable orientations of magnetization. 



expansion coefficient of the PZT stack (+1 ppm/K along 
and -3 ppm/K perpendicular to the piezoelectric stack). 
By applying a finite voltage to the piezoelectric during 
cooldown we are able to control the thermally-induced 
stress within Ae < 10"^ for T < 50 K. 

Application of ±200 volts on the piezoelectric results 
in a shift of the angle of magnetization switching by 
A(pH ~ 10°, which is comparable to the size of the 
hysteresis loop in TrAMR vs (pH scans at Vpzt — 0, 
Fig. [H^bjC). Applying a magnetic field of = 50 mT 
oriented at ipn — 62°, which corresponds to the mid- 
dle of the hysteresis loop, compensates for the thermally 
induced uniaxial strain anisotropy and restores the orig- 
inal degeneracy between [010] and [100] magnetization 
directions of the unstrained GaMnAs for Vpzt ~ 0. An 
additional strain is then applied by varying voltage on the 
piezoelectric. In Fig. [3] TrAMR is plotted as a function of 
measured strain Ae, the corresponding Vpzt are approxi- 
mately marked on the top axis (there is a small hysteresis 
in Ae vs Vpzt)- As additional compressive strain is ap- 
plied along [010], this direction becomes the easy axis of 
magnetization and magnetization aligns itself with [010]. 
As additional tensile strain is applied along [010] direc- 
tion, the [100] direction becomes the easy axis and po- 
larization switches by 90°. The switching occurs in a few 
steps, indicating a few-domain composition of our device. 
At Vpzt — the magnetization has two stable orienta- 
tions, M||[l00] and Af||[010], and the orientation can be 
switched by applying a negative or a positive voltage on 
the piezoelectric. Thus, the device performs as a bi-stable 



non-volatile magnetic memory with electrostatic control 
of the state. 

The center of the loop can be shifted by adjusting 
(fiH- e.g., a 1° change shifts the center of the loop by 
Ae « 3.5 • 10^^. As H increases, the size of the hysteresis 
loop decreases, and the hysteresis vanishes for H > 100 
mT; see inset in Fig.|3l At _ff < 40 mT the loop increases 
beyond the ±200 piezoelectric voltage span. In our ex- 
periments the magnetic field balances the residual strain 
due to anisotropic thermal expansion of the PZT. Alter- 
natively, intrinsic piezoelectric properties of GaAs can be 
utilized, in this case there will be no thermally induced 
strain and electrostatic switching of the magnetization 
direction can be realized without applying an external 
compensating magnetic field. Scaling of the piezoelectric 
element from 0.5 mm down to ^ 1 fim. will compensate 
for a small strain coefficient in GaAs ((^33* ~ 10 •df^^'^'^ at 
30 K), reduce operating voltage to a few volts and allow 
electrostatic control of individual memory cells. 

We model the strain along [100] and [010] direc- 
tions as an extra magnetostatic energy density term 
2e[ioo]K,sm^iipm + 45°) + 2e[oio]ife sin^((/3„ - 45°) = 
AeKg sin(2(/3m) + const. Then, for a single domain mag- 
net the free energy density can be written as 

+HM cos{(p.m - Vh) + Ae/^e sin{2(p,n) (2) 

omitting constant offset, where Ki, and are cubic, 
uniaxial and strain anisotropy constants; H is the applied 
in-plane magnetic field; and and Lpu are the angles 
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between [110] direction and magnetization and magnetic 
field respectively; see schematic in Fig. [TJ We assume 
that is the same for [100] and [010] directions. 

In equilibrium the magnetization orientation min- 
imizes the free energy, dE/difm — and d?E/dip'^ > 0. 
The TrAMR can be calculated from Eq.[l]for a given an- 
gle (fiH of the external field H. From the fits to the exper- 
imental TrAMR data we can extract the anisotropy con- 
stants Ki , Ku and ifg . The model captures all the essen- 
tial features of the data, and corresponding fits are shown 
in Fig.|3]for the strained and unstrained devices. For un- 
strained device anisotropy fields 2Ki/M — 40 mT and 
2Ku/M = 6 mT. These values are significantly smaller 
than the previously reported values for as-grown (not an- 
nealed) wafers^'-?'. For the sample attached to the piezo- 
electric the crystalline anisotropy field remains the same, 
but the uniaxial anisotropy increases to IK^jM = 50 
mT. The strain-induced anisotropy field /^eK^/M varies 
between 13 mT and 19 niT for different Vp^t between -200 
V and 200 V, the coefficient K^/M = 17 T. 



To illustrate the mechanism of magnetization switch- 
ing we plot magnetic energy density normalized by mag- 
netization E /M (Eq. [2]) as a function of in Fig. [S] 
At if = there are only two minima along the [100] axis 
due to the large uniaxial strain (see Fig. [5] (a, e)), caused 
by anisotropic thermal expansion coefficient of the piezo- 
electric. 

With external magnetic field ii = 50 mT applied at 
ifiH — 62° E/M has two minima: at fm = 32° and at 
123°, i.e., close to [010] and to [100] crystallographic di- 
rections. For the strain field Aeifg/M = 13 mT the 
global minimum is at = 32°, and in equilibrium the 
magnetization is oriented along the [010] direction. As 
the strain field increases to 19 mT the two minima switch, 
and ipm = 123° becomes the global minimum. It is inter- 
esting to note that there is always a small barrier between 
the two minima. Unless the barrier is an artifact of our 
model, the switching of magnetization should be either 
temperature activated or should occur via macroscopic 
quantum tunnelling. 
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